KEK-CP-072 
KEK Preprint 97-218 
November 1997 

K + — > 7r + 7T° Decay Amplitude with the Wilson Quark Action in 

Quenched Lattice QCD 



JLQCD Collaboration 
S. Aoki a *, M. Fukugita b , S. Hashimoto ct , N. Ishizuka a,d , Y. Iwasaki a,d , K. Kanaya a,d , 
Y. Kuramashi 6 , M. Okawa 6 , A. Ukawa a , T. Yoshie a,d 
3 Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
h Institute for Cosmic Ray Research, University of Tokyo, Tanashi, Tokyo 188, Japan 
'Computing Research Center, High Energy Accelerator Research Organization (KEK), Tsukuba, 

Ibaraki 305, Japan 

d Center for Computational Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
e Institute of Particle and Nuclear Studies, High Energy Accelerator Research 
Organization (KEK), Tsukuba, Ibaraki 305, Japan 
(February 7, 2008) 



Abstract 

We present a calculation for the K + — * tt + tt° decay amplitude using 
a quenched simulation of lattice QCD with the Wilson quark action at 
j3 = 6/<? 2 = 6.1. The decay amplitude is extracted from the ratio, the 
K — > 7T7T three-point function divided by either K and tt meson two-point 
functions or K meson two-point function and 1 = 2 irn four-point function; 
the two different methods yield consistent results. Finite size effects are exam- 
ined with calculations made on 24 3 x 64 and 32 3 x 64 lattices, and are shown 
that they are explained by one-loop effects of chiral perturbation theory. The 
lattice amplitude is converted to the continuum value by employing a one- 
loop calculation of chiral perturbation theory yielding a value in agreement 
with experiment if extrapolated to the chiral limit. We also report on the 
K meson B parameter Bk obtained from the K + — ► 7r + 7r° amplitude using 
chiral perturbation theory. 
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I. INTRODUCTION 



Despite the full understanding of the fundamental theory of weak interactions, the non- 
leptonic decay of hadrons still remains to be the least understood weak process, the most 
notable problems being AI = 1/2 rule and the calculation of e'/e. The predicament origi- 
nates from the difficulty of evaluating the hadronic matrix element of the product of currents. 
Much work has already been done to attack this problem using lattice QCD simulations |l] Q . 
but they have not yielded satisfactory results. 

Difficulties have proven to be especially severe for the the AI = 1/2 amplitudes |l] [|. 
From the computational point of view the problem lies in a calculation of the so-called eye 
diagram, which suffers from extremely large statistical fluctuations |l]-§,^J. Theoretically, 
this may be related to mixing of the dimension six weak operator responsible for the decay 
with operators of lower dimensions whose coefficients diverge linearly in the continuum limit. 
At a more fundamental level, there is the difficulty || that the K — > ttti 3-point function 
evaluated in Euclidean space-time does not yield information on the phase of the decay 
amplitude. 

Calculation of the AI = 3/2 process is known to be easier than that of the AI = 1/2 
process. For this case mixing of lower dimension operators is absent, and the so-called figure- 
eight diagrams which represent the AI = 3/2 amplitude, have clear signals in numerical 
simulations. Indeed lattice calculations have been reported by several groups 0,|J quite a 
long time ago. The problem, however, was that the results turned out to be inconsistent 
with experiment: lattice calculations have given the amplitude roughly a factor of two larger 
than experiment. 

Two potential origins are suspected to give this discrepancy. One is an issue in the 
matching of the lattice and continuum operators. Early studies employed the factor \^2k for 
the quark wave function normalization and the bare lattice coupling constant for estimating 
the renormalization factor of the four-quark weak operator. It is by now well known that 
the KLM factor yl — 3k/4k c and tadpole- improved perturbation theory || are more 
adequate for the operator matching. 

Another problem concerns the use of chiral perturbation theory (CHPT) to convert lattice 
results into the physical amplitude. Only the tree-level formula was known and used in the 
previous work. The meson mass dependence of lattice calculations appeared consistent with 
the prediction of the tree-level formula, allowing, however, for large statistical errors. It was 
probably necessary to use the formula including higher order CHPT effects, but its necessity 
was not manifest. An interesting development in this connection is a recent calculation of 
one-loop corrections to the K + — > 7r + 7r° amplitude in CHPT by Golterman and Leung j9]. 
Applying their results to the old data obtained by Bernard and Soni ||, they found that 
one-loop effects decrease the physical amplitude by about 30%. 

With the hope to improve the problems posed here, we have carried out a high statistics 
simulation of the K + — > 7r + 7r° amplitude in quenched lattice QCD, incorporating various 
theoretical and technical developments made in recent years. In particular, we discuss in 
detail how one-loop corrections of CHPT affect physical predictions for the decay amplitude 
from lattice QCD simulations. We also report on the K meson B parameter Bk obtained 
from the K + — > 7r + 7r° amplitude using CHPT. 

This paper is organized as follows. After a brief description of simulation parameters in 
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Sec. m, we explain our method for extracting the decay amplitude in Sec. |TJ Our results 
for the K + — > 7r + 7r° amplitude are presented in Sec. [IV| with discussion made on one-loop 
effects of CHPT. Results for Bx are given in Sec. M. Sec. [VT| summarizes our conclusions. 



II. SIMULATION PARAMETERS 

Our simulation is carried out in quenched lattice QCD employing the standard plaquette 
action for gluons at (3 — 6.1 and the Wilson action for quarks. We take up, down and strange 
quarks to be degenerate, and make measurements at four values of the common hopping 
parameter, k = 0.1520, 0.1530, 0.1540, and 0.1543, which correspond to M^/M p = 0.797, 
0.734, 0.586 and 0.515. In order to examine finite-size effects, simulations are carried out 
for two lattice sizes, 120 configurations on 24 3 x 64 and 65 configurations on 32 3 x 64. 
Gluon configurations are separated by 2000 pseudo heat bath sweeps. Quark propagators 
are solved with the Dirichlet boundary condition imposed in the time direction and the 
periodic boundary condition in the space directions. 

We adopt 1/a = 2.67(10)GeV for the physical scale of lattice spacing estimated from the 
p meson mass, and k c = 0.15499(2) for the critical hopping parameter, which were obtained 
in our previous study [P~D| . Our calculations are carried out on the Fujitsu VPP500/80 
supercomputer at KEK. 



III. METHODS 
A. Extraction of decay amplitude 

Let us consider the four-quark operator defined by 

Q+ = 2 • [ - 7s)d WTAi(! - l5)u + sj^l - j 5 )u uj^l - 75)^ ] , (1) 

which is relevant to AI = 3/2 two-pion decay of the K meson. We first discuss our method 
for extracting the lattice matrix element of the operator Q + , deferring the question of match- 



ing the lattice and continuum operators to Sec. II 1 B 



We extract the decay amplitude from the 4-point correlation function 

G Q (t + ,t ;t;t K ) = (0\W + (t+)W (t )Q + {t)W K (t K )\0) . (2) 

In order to enhance signals we construct wall sources (denoted by W) for all external mesons, 
and fix gauge configurations to the Coulomb gauge. The wall sources Wk, W+, and Wo for 
K + , 7r + , and 7r° are placed at the time slices tx, t + and to such that £r- <C t <S t+, to- All 
mesons are at rest, and the 4-quark operator Q + is projected to zero spatial momentum. 

In our calculation for temporal lattice size T = 64, we place the K meson at tx = 4. The 
two 7r mesons are placed at different time slices, t + = 59 and t = 60 to avoid contaminations 
from Fierz- rearranged terms in the two-pion state that would occur for the choice t + = to. 

The correlation function Gq behaves for tx <C t t + ~ to as 
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Gc}(t+, t ; t; tx) 

= (0|^ + (0)^ (t - t + )|7r+7r°)-^(7r + 7r°|g + (0)|^ + >-L(^ + |W^(0)|0) 



K 



where denotes the normalization factor of the K meson state, |7r + 7r°) represents the 
1 = 2 two-pion state with a mass M n7T and a state normalization factor N n7T . 

In order to remove the normalization factors in Gq we calculate the product of the meson 
2-point functions given by 

G w (t + ,t ;t;t K ) = (0|iy (to)vr°(t)|0)(0|iy + (t + )7r + (t)|0)(0|^ + (t)^(^)|0) . (4) 

Defining a ratio R\y = Gq/Gw, we find 

Rw(U, to; t; t K ) = S W ■ ^MnJf +) ' e(t ~ t+)A ( 5 ) 

where 

A = - 2M 7r (6) 
is a mass shift due to a finite spatial lattice size, and Sw is defined by 

« = jg. (0|^ + (0)^o(to-t + )|7r+7r°) 
w (0|W + (0)|7r+)(0|W (to-t + )K°> lJ 

where to — 1 + = 1 in our calculation. The value of Sw should converge to unity for infinite 
volume. 

In Fig. D we plot (7r|7r|0) 3 ■ Rw at k = 0.1530 as a function of time t of the weak operator, 
where we calculate (7r|7r|0) from the pion 2-point function for point source and point sink. We 
observe a clear non- vanishing slope, which means the mass shift A being positive. Numerical 
values of A and the decay amplitude (tt + tt°\Q + \K + ) obtained by a single exponential fit for 
the time range t = 18 — 46 are tabulated in Table |. Here we assume Sw = 1? whose 
justification will be discussed below. 



According to Luscher's formula ||11|| , the finite-size mass shift of the two-pion state is 
written 

A = M ra -2M, = -^ + (i -) (8) 

where a® is the s-wave scattering length and L is the spatial size. This formula was previously 
employed to calculate the s-wave im scattering length in quenched lattice QCD fl2|-|i"4T| . It 
was found that lattice calculations give ao in good agreement with the prediction of current 
algebra. Using the current algebra formula ap =2 = M n /(8irF^) with F n = 132MeV and 
l/a = 2.67(10)GeV, we obtain aA = 0.015 for L = 24 and 0.006 for L = 32. Considering 
uncertainties arising from terms of 0(L~ A ) and the difference between the physical and 
measured values of F n , we regard this estimate being consistent with the measured a A (see 
Table H). 
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As an alternative method we may remove the normalization factors of the 4-point function 
Gq with 

G P (t + ,t Q ;t;t K ) = (0|^ + (t + )^ (t )7r + (t)7r°(t)|0)(0|K + (t)^(^)|0) . (9) 

The ratio Rp = Gq/Gp is independent of t and it does not depend on the wall sources for 
t K < t < t + ~ t , 

Rp(t+, t ; t; t K ) = Sp 1 ■ - — - J^jo — - (10) 

{'n-\'n-\U)- > 

where 

(tt+ttV^O) 

^ P = (vrlvrlO)^ (U) 

which should become unity for infinite spatial lattice. 

The dependence of Rp on the time t of the weak operator is shown in Fig. |2] for k = 0.153, 
the same hopping parameter as in Fig. [1] for R w . As expected, a clear plateau is seen for 
t « 20 — 40, where effects of excited states near the lattice boundaries already disappear. 

In Table | we list (7T +, k \Q + \K + ) obtained by fitting Rp to a constant over t = 22 — 42 
assuming Sp — 1. The results from the two methods show good mutual agreement, well 
within the statistical error of 10 — 15%. We note that statistical errors for Rp are smaller, 
and therefore adopt the matrix elements from Rp to obtain the physical decay amplitude 
below. 

We still have to justify the assumption Sw = Sp = 1 used above. This is not a priori 
obvious, especially for Sw, since wall sources are uniformly extended across the spatial 
lattice, although a good agreement of (tc + tc \Q + \K + ) from Rw and Rp implies Sw • Sp close 
to unity. For Sp chiral perturbation theory predicts a finite-size correction of the form || , 

M 2 

Sp = 1 + 24F*{M n aL)* ' (12) 

This formula indicates that the deviation of Sp from unity would be less than 1% in our 
simulation. Hence we conclude Sw ~ 1. 



B. Operator Matching 

For the quark field normalization we employ the KLM factor 



/ 3k 

^continuum ^lattice (13) 

V 4ac c 

Due to CPS symmetry the weak operator Q + defined in (]!]) does not mix with other op- 
erators [0]. With the tadpole improvement (with the factor uq = 1/8k c ) the multiplicative 
renormalization factor for Q + , which relates the lattice operator to the continuum one at a 
scale /i, is given by |il||T6| 
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Z(ai) = 1 + ^^[-41og(//a/7r) - 21.140] (14) 

where the naive dimensional regularization (NDR) is taken with the modified minimum 
subtraction scheme (MS) in the continuum. 

We employ the MS coupling constant estimated as follows. First we obtain g v by [f| 

- logP = i^(3.41/o){l - (1.19 + 0.017JV/) g y^ /a) + 0(4)} (15) 

with P the average plaquette. Next we calculate Ay from g v using 

3.41/a 2 1 3i flix/do / 1fi x 

log( — > = ^ V° g TT7WA (16) 

where O = 11 - 2N f /3, Q x = 102 - 38A^/3, and x = # y (3.41/a)/(47r) 2 . A perturbative 
relation A^ = 0.6252 • Ay then yields Aj^r, with which we can calculate <7fjg(/f) & t an Y 
scale [A. In the present calculation we find Aj^g = 293(ll)MeV with P = 0.605 and 1/a — 
2.67(10)GeV at 3 = 6.1. 

Let A2 be the physical amplitude for AI = 3/2 K — > tttt decay. Experimentally, 

y/| • ReA 2 • • V* s V ud }- 1 = 10.4 x lO" 3 GeV 3 . (17) 
The relation of the decay amplitude to the matrix element of Q + is 

y| ■ ReA 2 ■ ■ CKJ- 1 = cf'V) ■ (n+AQ ( + Nf \»)\K + ) . (18) 

On the right-hand side C^ f \fx) and Q+ (/J,) are the Wilson coefficient function and the 
renormalized weak operator at a scale fi with superscript Nf the number of quark flavors 
appropriate for the scale \i. We choose /1 = 2GeV to estimate the physical amplitude, and 
hence Nf = 4. 

In our calculation, matching of the lattice operator Q l ^ tUce to the continuum operator 
(2GeV) is not straightforward since the simulation is carried out in quenched QCD 
(Nf = 0). To treat this problem we proceed in the following way. We first match the lattice 
operator to the continuum operator for Nf = at a scale q* using the renormalization 
factor Z(q*) in flnj): Q { +\q*) = Z(q*)Q l f tice . The operator Q { +\n) at any scale fj, can then 
be obtained by renormalization group evolution in the continuum: 

= U {0) (fi,q*)Z(q*)Q l f ice (19) 

where U^ Nf \jj, fx') is the two-loop renormalization group running factor from scale // to /i 
and it is given by 
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Here 70 = 4 and 71 = — 7 + 4iVy/9 are the one- and two- loop anomalous dimensions for Q + 
III- 

In the spirit of tadpole improvement, the matching point q* from the lattice to the 
continuum operator should be chosen to minimize higher order contributions in the renor- 
malization factor Z(q*). Since an estimate of this value is not available, however, we take 
q* = 1/a or ir/a and investigate the q* dependence of the decay amplitude. 

We still need to relate the operator of the Nf = theory to the operator of the 
Nf = 4 theory. Whether such a matching is possible is a problem generally encountered in 
quenched QCD calculations of weak matrix elements. As a working hypothesis, we assume 
that there is a scale k*, typical of the K + — > process, at which the Nf = operator 

matches with the Nf = 4 operator, 

U^(k*,q*)Qf\q*) = Q? ) (F) . (21) 
We then estimate the decay amplitude for the Nf = 4 theory by 

= cf (/x) ■ U^(fi,k*) ■ (n + n°\Q^(k*)\K + ) 

= CfV) • D(fi, k*, q*) ■ (ir+7r°\Q l f ice \K + ) (22) 

where 

D(fx,k*,q*) = U {4 \^k*)U {0) (k*,q*)Z{q*) . (23) 

For the renormalization group evolution in the continuum we follow Buchalla et al. |18| . In 
particular we use their cl 4) (2GeV) = 0.859 with AgL = 215MeV. 

The value of the matching scale k* is not known. The variation of D(fi,k*,q*) with 
respect to the scale k*, however, arises from the difference of the A parameter and the 
anomalous dimension of Q + for Nf = and 4, and so it is expected to be small. The values 
of D(fi, k*,q*) for several k* are tabulated for q* = 1/a and n/a in Table [0]. We observe that 
the dependence on k* is indeed very small, and we set k* = lGeV in the following analysis. 

Let us note that the difference of D(fi,k*,q*) for q* = 1/a and ir/a is about 10%. This 
is the largest systematic error in our operator matching procedure beside the assumption of 
the matching scale k*, and it is comparable to our statistical errors. 



IV. RESULTS FOR THE K+ -> vr+vr AMPLITUDE 
A. Decay amplitude with tree-level CHPT 

As in the previous work 0Qj we take degenerate strange and up-down quarks, and 
assume all external mesons at rest. The amplitude obtained with this kinematics is clearly 
unphysical, having an energy injection at the weak operator. In order to relate the lattice 
result to the physical amplitude information is needed on the dependence of the amplitude 
on the K and n masses away from the physical point. 

Earlier calculations have used chiral perturbation theory (CHPT) at tree level for this 
purpose. The operator Q + is decomposed under chiral SU(3) L into terms belonging to 
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[8, AI = 1/2], [27, AI = 1/2] and [27, AI = 3/2]. The [27, AI = 3/2] part of Q+, which 
contributes to fT + — > 7r + 7r°, is given by 

I ■ Q 4 = ~ ■ [2 ■ Q + - s 7ll (l - l5 )d d 7ll (l - 75 )d] . (24) 
In general the 27 operator in QCD can be described by operators in CHPT as, 

2 Q 7 CD = a 27 ■ R% • (E^St) ife (S^St) iZ (25) 

where, for Q 4 , the non- vanishing components of the tensor are R^\ = R\% = R^ = 
R%\ = | and R%% = R%3 — ~ \ > an d the pseudoscalar meson field is given by 

S = e i7r// (26) 

for the full theory, or 

S = e^/'e*" 7 ^ (27) 

for the quenched theory. At tree level of CHPT one obtains the formula connecting the 
physical amplitude and that calculated on the lattice 0] : 

(,+AQ+\K + ) phys = • (J) • (^) 3 ■ (^[Q+I^We (28) 

where mx = 497MeV and = 136MeV are physical masses, and M n is the degenerate 
K and n masses on the lattice. We emphasize that the constant a 2 7 and the tree-level 
decay constant / may take different values in the full and quenched theories. We denote the 
constants in quenched theory with superscript q. 

In Fig. ^| we compare the decay amplitude C + ■ (ti + 'k \Q + \K + ) of the previous work at 
(3 = 5.7 and 6.0 [0] with ours at f3 — 6.1. Here, as a working hypothesis, we set a 27 and / to 
be equal in full and quenched theories. For the sake of comparison, our data are analyzed in 
a manner parallel to that in Ref. || as much as possible, i.e., employing the traditional v2k 
normalization for quark fields, no tadpole improvement in the renormalization factor, and 
applying the tree-level relation (|28|). The matching factor ( p3| ) with q* = nja is applied in 
our results for consistency. Since the normalization adopted in Ref. P| for comparison with 
experiment differs from ours, we plot the results divided by the experimental value. In view 
of various differences in the simulation parameters and details of analysis procedures, the 
values from the two studies are taken to be consistent, both being larger than experiment 
roughly by a factor of two. 

Let us note that our results, which attain errors of about 10%, show a clear dependence 
on the lattice meson mass M n . The presence of finite-size effects is also evident, exhibiting 
a decrease between 24 3 and 32 3 spatial sizes. We observe that these features were present 
in the previous simulations when examined in the light of our results, but they were not 
evident at the time because of large statistical errors of 20 — 30%. 

In Fig. | we show how the use of the KLM normalization affects the meson mass depen- 
dence of the decay amplitude (plotted with filled symbols) as compared to the conventional 
normalization ^2k/2k c /2 (open symbols). While the amplitudes for small M w change only 
slightly, those for larger M n increase by about 20%, which is beyond the statistical error by 
a factor of two. A significant meson mass dependence and finite-size effects observed in our 
data show that tree-level CHPT is inadequate to extract the physical amplitude from lattice 
calculations. 
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B. Decay amplitude with one-loop CHPT 



Recently Golterman and Leung have carried out a one-loop calculation of CHPT for 
the decay amplitude in full and quenched QCD for degenerate and non-degenerate K and ir 
mesons . Their formula also includes finite-size correction terms. Combining with the one- 



loop formula calculated for the physical point ||19|| , we analyze how our results are changed 
by one- loop effects of CHPT. 

Let us denote by (n + 'K \Q + \K + ) p hys the physical amplitude in the full theory with non- 
degenerate K and 7r mesons of mass rrix and m n , and by (n + n°\Q +\K + ) i attice the amplitude 



in the quenched theory with degenerate K and tt mesons of mass M n 
man and Leung, 



According to Goiter- 



(7i + n°\Q + \K + ) 



m 



phys 



K 



2MI 



I /«27\ (U 

" WW 



■Y-(7r + 7i°\Q + \K + ) 



lattice 



(29) 



where a and / are defined in (2~5-E7|), and the factor Y is given by 



Y 



1 + 



U + d 



1 + 



Ml 



(47rF^) 2 



-31og(ff) + F{M n aL) + d q 



(30) 



The numerator of Y represents the one-loop effect in the full theory, and the denominator is 
the corresponding effect in the quenched theory. The dimensionless constants d and d q are 
the contact term coefficients arising from the 0(p 4 ) terms of the chiral Lagrangian. f n and 
F w are the one-loop corrected decay constants in the full and quenched theories, which differ 
from the tree-level values / and f q . In the numerator of Y, U is a complicated function of 
physical K and 7r masses, the decay constant /„. and fx, and the cutoff of CHPT for the 
full theory A coni , and a numerical approximation is 



U — A 



B ■ log( 



A 



cont 



(31) 



where A = -104.73 and B = -29.57 for m n = 136MeV, m K = 497MeV, f w = 132MeV, and 
f K = 160MeV. In the denominator of ©, A q is the cutoff of CHPT for quenched QCD, 
and F(M 7r aL) represents finite-size corrections for a spatial size L which takes the form 



F(M v aL) 



17.827 



+ 



12vr 2 



M n aL (M^aLf 



(32) 



We set «27 and / to be equal in the quenched and full theories as in the analysis with the 
tree-level CHPT. We initially ignore the effects of 0(p A ) terms of the chiral Lagrangian d 
and d q . We leave A cont and A q to be different, however, and examine the dependence of the 
results on these cutoffs. 

In Fig. H we plot the one- loop corrected decay amplitude for A q = 770MeV and lGeV 
for the choice A cont = 770MeV. We set f w = F n = 132Mev in Q3IJ), and the finite-size 
corrections F(M n aL) are taken into account. The results of a similar analysis for the choice 

jycont = 1Ge y are p l ottec i j n Fig _ g 

An important feature observed in Figs. [5] and |6] is that the size dependence seen with the 
tree-level analysis in Fig. 01 is removed after finite-size corrections at the one-loop level. At 
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the same time, the amplitude decreases by 30 — 40% over the range of meson mass covered 
in our simulation. 

Another noteworthy feature in Figs. [| and |6] is that a sizable lattice meson mass depen- 
dence still remains in the amplitude, and that the magnitude of the slope depends sensitively 
on the choice of A q . This feature can be understood as arising from the 0(p 4 ) coupling con- 
stants in the quenched theory, i.e., d q in fl3~0D , which was ignored above. If we denote our 
present results by (n + TT \Q + \K + ) ours we find from (]3"0|) that 



11 M i d 

1 * ( A -~ TP \2 U 



7T 



-AQ + \K + )ours = (7r + AQ + \K + ) phys (33) 



1 J m * A 

1 + (M,P a 



showing the presence of a term linear in M 2 . We note furthermore that d q actually depends 
on A 9 : d q = d q (A q ). Since the total 0(p A ) correction in the denominator of (|30| ) should be 
independent of the cutoff A 9 , d q for different values of A q varies according to 

^(A 9 )=d,(A' 9 )-31og(^) 2 . (34) 

To compare these relations with our results, we fit ('K + 'n' \Q+\K + ) ours as a function 
of M 2 to the form ([33]). Employing data for which the value of M n does not exceed 
the cutoff, we find ^(lGeV)/(4vr) 2 w 0.015 and rf g (770MeV)/(47r) 2 w 0.025. The differ- 
ence d 9 (lGeV)/(47r) 2 — d g (770MeV)/(47r) 2 « —0.01 is in good agreement with the value 
-3 log (lGeV/770MeV) 2 /(4vr) 2 = -0.0099 expected from These results show that the 
uncertainties associated with d q can be removed by a chiral extrapolation of our amplitude 
to the chiral limit M n = 0. 

A further consequence of (^) is that a correction due to the 0(p 4 ) term d in the full 
theory remains even after taking the limit M n — > in our results. In order to examine the 
magnitude of this uncertainty, we use an estimate d(A cont ) / (An) 2 = 0.003(14) at A cont = m v 
from a phenomenological analysis |19[ . In view of the formula 



/ A cont \ ^ 

d(A c ° nt ) = d{A' c ° nt ) - 29.57 • log ( —J (35) 



obtained from (H) and (|3T|), this leads to a value ci(770MeV)/(47r) 2 w -0.12 and 
d(lGeV) / (Ait) 2 w —0.22. These values imply that the physical decay amplitude is 10% 
lower than our results for A cont = 770MeV, and 20% for A cont = lGeV. This provides an 
explanation of a discrepancy of about 10% observed in Figs. [| and || between the values of 
(n + n°\Q + \K + ) ours calculated with A cont = 770MeV and lGeV. Let us add a remark that 
the values of d estimated above for A cont = 770MeV and lGeV is an order of magnitude 
larger compared to those of d q for the quenched theory. 

We find from this analysis that including the correction due to the 0(p A ) coupling con- 
stants is possible if an accurate value of d is known from phenomenological studies. Since 



this is not yet the case |19|, we shall not pursue this point further here, leaving the correction 
as a source of uncertainty in our final results. 

The amplitude obtained from (vr + 7r°|Q + |7^ + ) ou , rs by a chiral extrapolation to the limit 



M n = is listed in Table 111 for several choices of the cutoff and the operator matching 
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point q*. In the results in Table [TJ, the systematic error due to the matching scale q* is 
about 10%. Statistical errors are larger (about 20%), mainly due to a linear extrapolation to 
the chiral limit. Within these uncertainties and that of 10-20% due to the d term discussed 
above, the values in Table ITXTI are consistent with the experiment 10.4 x 10 _3 GeV 3 . 



V. B K FROM THE K+ -> vr+vr AMPLITUDE 

The AS = 2 four-quark operator defined by 

A s=2 = S7^(l - 7 5 )d S7 M (1 - 7 5 )c? (36) 

belongs to the same 27 representation as the operator Q 4 which is the [27, AI = 3/2] part of 
Q + . As a consequence one can obtain the K meson B parameter Bk from the K + — > 7r + 7r° 
amplitude using CHPT. The non-vanishing component of the tensor i?^ in (^) for this 
operator is given by = 1. 

The one-loop relation in CHPT for quenched QCD for the unphysical degenerate case 
has been obtained by Golterman and Leung M, 

_J_ 3-(7r + 7r°|Q + |A- + ) 



with the one-loop correction R given by 

R 



M 2 r M -i 

f 31og O 2 + F(M.aL) . (38) 



(4^ 



;2 



Here F# and F„. denote the K and 7r meson decay constants in the quenched theory, and 
the other notations are the same as those in (|30H32|). 

Our procedure for calculating Bk from fl37|) is essentially the same as for the K + — > 7r + 7r° 
amplitude including the operator matching procedure, although the coefficient function C + 
is absent in the present case. In Figs. and || we plot _B^(2GeV) obtained from the 
K + — > 7T + 7r° decay amplitude with tree and one-loop CHPT. We set = 160MeV and 
F n = 132MeV in Q57[ ) and fl3"8p. The one- loop CHPT effect and the cutoff dependence for 
a small M% region are small compared with those for the decay amplitude. At the physical 
K meson mass Ml = 0.246GeV 2 B K takes almost the same value for different choices of 



A q and the lattice size. In Table |TV| the average of the two data points with the smallest 
M w is tabulated. Our results, B K = 0.581(56) — 0.663(67) are consistent with the JLQCD 
value £>x(2GeV) = 0.68(11) POfl obtained at the same coupling constant (3 = 6.1 through a 



calculation of the K° — K matrix element of the AS = 2 operator employing chiral Ward 
identities for determining the mixing coefficients. 

A direct calculation of Bk with the Wilson quark action has the complication that 
the operator mixing problem of the AS = 2 operator has to be solved non-perturbatively, 
which causes large statistical errors. In contrast, the Q + operator does not mix with other 



operators as mentioned in section [IIIB| . Therefore, statistical errors of Bk obtained from 
the K + — > 7r + 7T° amplitude is smaller. Theoretical uncertainties associated with the use of 
CHPT, however, are large in this approach that offsets the advantage of the present method. 
In any case, our calculation, albeit with a significant error, provides an independent check 
for Bk for the Wilson quark action obtained with the chiral Ward identity procedure, and 
also supports the validity of CHPT. 
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VI. CONCLUSIONS 



In this article we have reported results of a study of the K + — > 7r + 7r° decay amplitude in 
quenched lattice QCD. With a set of high statistics simulations we have found that the results 
show sizable finite-size effects, which, however, are consistent with those predicted by a recent 
one-loop calculation of CHPT. We have furthermore seen that a meson mass dependence 
which remains after inclusion of the one-loop corrections of CHPT in the prediction for 
the decay amplitude is due to effects of the 0(p 4 ) contact terms in the quenched theory. 
Making an extrapolation to the chiral limit to remove these effects, we have found 8.9(1.7) x 
1CT 3 — 11.4(1.5) x lCT 3 GeV 3 for the physical value of the decay amplitude, depending on 
the choice of the cutoff parameter of CHPT. These values are consistent with experiment 
(10.4 x l(T 3 GeV 3 ). 

The present result may be compared to those of the previous studies p,|J which gave 
decay amplitudes roughly twice larger than experiment. Our smaller value originates from 
the two effects, one-loop corrections as also noted by Golterman and Leung in their reanalysis 
of the old results, and a decrease of the amplitude toward smaller values of M n . 

As a further application of the one-loop formula, we have calculated the Bk parameter, 
and found that it is consistent with a recent direct calculation for K — K mixing. 

The encouraging results we have obtained, however, should be taken with several reser- 
vations. The value of the K + — ► 7r + 7r° decay amplitude estimated in the chiral limit suffers 
from uncertainties of 10-20% due to the 0(p 4 ) contact terms of the full theory, because the 
phenomenological estimate available is not very accurate. A sizable finite-size correction of 
30—40%, while consistent with the one-loop prediction of CHPT, raises the question whether 
ignoring higher order corrections can be justified. Furthermore, various constants of CHPT, 
in particular the coefficient «27, may differ between the quenched and full theories, and we 
have no way of estimating or correcting the difference. Reliability of CHPT for calculating 
unphysical amplitudes could also be an issue. Reducing these sources of uncertainties, es- 
pecially those related to quenching and better controlling finite-size effects would require a 
difficult task of carrying out simulations in full QCD on a physically large lattice. 
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FIG. 1. (7r|7r|0) 3 • Rw{t + , to; t; tx) at k = 0.153. Open and filled circles refer to data for 24 3 
and 32 3 lattices. 



14 



<7C|7C|0> 3 X FL(t) 



k = 0.1530 (IVI^/M = 0.734) 
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FIG. 2. (7r|7r|0) 3 • Rp(t + ,t ;t;t K ) at k = 0.153. Open and filled circles refer to data for 24 3 
and 32 3 lattices. 
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FIG. 3. Comparison of our results for C+-(n + Tr°\Q+\K + ) normalized by the experimental value 
obtained with tree-level CHPT relation (^) for q* = ir/a at (3 = 6.1 with those of previous work 
0] at (3 = 6.0 and 5.7. Results are plotted as a function of lattice meson mass M%. Traditional 
\^2k normalization is employed for quark fields and tadpole-improvement is not applied in the 
renormalization factor. 
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FIG. 4. Decay amplitude C+ • (it + it \Q+\K + ) for q* = 1/a as a function of lattice meson mass 
M% for tree-level CHPT. Circles and squared refer to decay amplitude obtained on a 24 3 and 32 3 
lattice. Open symbols correspond to the traditional ^2k/2k c /2 normalization factor of Wilson 
quark fields, while filled symbols are for the KLM normalization. 
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FIG. 5. Decay amplitude C + ■ (ir + ir°\Q + \K + ) for q* = 1/a obtained with one-loop CHPT for 
A cont = 770MeV plotted as a function of M%. Circles and squares refer to data for 24 3 and 32 3 
spatial sizes. Open symbols are for A q = 770MeV and filled symbols for A q = lGeV. 
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FIG. 6. Decay amplitude C+ ■ (7r + vr |Q + |K + ) for q* = 1/a ELS cl function of obtained with 
one-loop CHPT for A con * = lGeV. Meaning of symbols are the same as in Fig. ||. 
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FIG. 7. -Bx(2GeV) for q* = 1/a obtained from K + — > 7r + 7r° decay amplitude as a function of 
M% obtained with tree level CHPT. Circles and squares refer to data for 24 3 and 32 3 spatial sizes. 
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FIG. 8. Bk (2GeV) for q* = 1/a obtained from K + — > 7r + 7r° decay amplitude as a function 
of M% obtained with one-loop CHPT for A q = 770MeV and lGeV. Circles and squares refer to 
data for 24 3 and 32 3 spatial sizes. Open symbols are for A q = 770MeV and filled symbols for 
AQ = iGeV. 
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TABLES 



K 






aA 

(io- 3 ) 




(ir+Tr°\Q + \K+) 
from Ry/ from Rp 


L = 24 

0.1520 

0.1530 

0.1540 

0.1543 


0.3440(14) 
0.2776(17) 
1967(19") 
0.1653(21) 




7.9(1.7) 
8.6(2.0) 
9 3(2 9) 
8.6(3.6) 




0.261(20) 
0.151(14) 
0.0617(81) 
0.0382(58) 


0.271(19) 
0.160(13) 
0.0680(69) 
0.0434(49) 


L = 32 
0.1520 
0.1530 
0.1540 
0.1543 


0.3459(10) 
0.2784(11) 
0.1914(13) 
0.1651(15) 




3.6(1.4) 
4.2(1.5) 
5.7(2.1) 
7.1(2.8) 




0.229(17) 
0.132(13) 
0.0565(71) 
0.0383(55) 


0.234(16) 
0.135(11) 
0.0573(55) 
0.0380(37) 


TABLE I. 

assume S\y = 
Rw and by a 


The mass shift A = M wn - 2M n and (7r+7r°|Q 
Sp = 1. These value are obtained by a single ( 
constant fit over t = 22 — 42 for Rp. 


4- \K + ) from Rw and Rp. Here we 
;xponential fit over t = 18 — 46 for 


fc*(GeV) 






q* = l/a 






q* = ir/a 


0.700 
1.000 
1.500 
2.000 






0.759038 
0.761556 
0.765198 
0.768126 






0.830913 
0.833670 
0.837657 
0.840863 


TABLE II. Values of D(2GeV, k* , q*) for A^ 


= 215MeV and A^L = 293MeV. 

MS 


^cont 

(GeV) 


A q 

(GeV) 


q* = l/a 


C+ ■ (7T+ 

24 3 

q* = 


AQ+ 

ir/a 


| J FC+)(xlO- 3 GeV 3 ) 

32 3 

q* = l/a 


q* = n/a 


0.77 
0.77 
1.0 
1.0 


0.77 
1.0 
0.77 
1.0 


9.3(1.9) 
9.4(1.3) 
10.3(2.1) 
10.4(1.4) 


10.2(2.1) 
10.3(1.4) 
11.3(2.3) 
11.4(1.5) 


8.9(1.7) 
8.8(1.1) 
9.8(1.9) 
9.7(1.2) 


9.7(1.9) 
9.6(1.2) 
10.7(2.1) 
10.6(1.3) 



TABLE III. Results of linear extrapolation of C + ■ (7r + vr°|Q + | K+) to Ml = 0. For 
A q = 770MeV fits are made with three points with smaller M w as M w of the forth point ex- 
ceeds the cutoff. Statistical and extrapolation errors are combined. The experimental value is 
10.4 x 10~ 3 GeV 3 . 
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tree 


A« = 0.77GeV 


Ai = lGeV 


24 3 


0.728(78) 


0.587(64) 


0.659(71) 


32 3 


0.627(63) 


0.581(58) 


0.663(67) 



TABLE IV. B K (2GeV) at physical K meson mass = 496MeV obtained from the 
K + — > 7r + 7r amplitude. The row "tree" refers to result with lowest CHPT and others are ob- 
tained by onedoop CHPT for M = 0.77GeV and lGeV. 
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